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ABSTRACT 


The  starting  process  of  supersonic  flow  in  a nozzle 
followed  by  constant  area,  straight  duct  diffusers  was 
studied  in  a Ludwieg  tube  equipped  with  an  upstream 
diaphragm.  The  time  required  to  establish  steady  flow 
from  rest  at  the  exit  of  the  supersonic  nozzle  was 
determined  from  shadowgraph  pictures  of  the  flow  field, 
and  from  static  pressure  measurements.  The  initial 
pressure  ratios  across  the  diaphragm  to  start  the  super- 
sonic flow  at  the  nozzle  exit  are  also  reported.  They 
are  compared  with  the  minimum  pressure  ratio  required  to 
sustain  the  flow  at  the  nozzle  exit  during  the  first  flow 
cycle  of  the  Ludwieg  tube  operation. 

The  flow  Mach  number  is  M=2.8  and  nozzle  exit  area 

. . 2 
and  diffuser  cross  section  are  (2x2)  in  . Various  lengths 

of  diffusers  were  used  with  L/D=0,  1.5,  3,  4.5,  6.  Within 

5 6 

the  range  (2.8x10  to  1.7x10  ) of  Reynolds  numbers,  Re^,  as 
based  on  the  nozzle  exit  hydraulic  diameter,  the  experimental 
results  were  unchanged.  Starting  times  at  the  nozzle  exit 
with  boundary  layer  suction  to  remove  2%  of  the  total  mass 
flow  rate  applied  through  a slit  at  the  nozzle  exit  were 
16-33%  lower  than  the  ones  obtained  without  suction,  for  a 
3 inch  long  diffuser.  For  the  same  diffuser,  the  limiting 
initial  pressure  ratio  across  the  diaphragm  that  allows  the 
smooth  start  of  the  supersonic  flow  at  the  nozzle  exit , was 
unchanged  by  suction. 
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NOMENCLATURE 


a speed  ?f  sound 

. 2 

A area,  nozzle  exit  area  =(2x2)  in 

D hydraulic  diameter  at  nozzle  exit  (4  times  area 

divided  by  perimeter),  =2  inches 
h*  height  of  the  nozzle  at  the  throat,  =0.472  inches 

Z characteristic  length  of  the  nozzle,  defined  as 

/ R*h*  = 1.89  inches 
L diffuser  length 

M Mach  number,  = v/a 

M shock  Mach  number 

_s 

M average  shock  Mach  number 

s 

p static  pressure 

R*  radius  of  curvature  of  the  diverging  section  of 

the  nozzle  contour  at  the  throat,  =7.64  inches 
Re^  Reynolds  number  based  on  hydraulic  diameter  and 

test  section  free  stream  parameters,  =vpD/p 
t time 

t time  interval  between  the  arrival  of  the  shock  at 

c 

the  nozzle  throat  and  the  occurrence  of  the  spark 

t total  starting  time,  =t  ,+t  „ 

s ^ s 1 s 2 

t dimensionless  starting  time,  =1  /t. 

S SI 

t ^ time  elapsed  from  the  diaphragm  rupture  until  shock 

wave  arrives  at  a specific  location 

t 0 time  between  the  arrival  of  the  shock  and  the 

s 2 

attainment  of  the  steady  free  stream  pressure 
T temperature 

v flow  speed 

x axial  distance,  nozzle  exit  x=0 

p density 


IV 


y dynamic  viscosity 

ij.  characteristic  flow  time,  defined  by  £/a*  = 0.153  ms 


Subscripts 


1 

4 

3 

o 

n 

d 


downstream  initial  conditions  before  diaphragm  rupture 
upstream  initial  conditions  before  diaphragm  rupture 
conditions  upstream  of  nozzle  after  the 
expansion  fan  has  passed  into  the  supply  tube 
nozzle  supply  conditions 
conditions  at  the  nozzle  end 
conditions  at  the  diffuser  end 


during 

first 

cycle 

of  the 

Ludwieg 

tube 

flow 
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I.  INTRODUCTION 

With  the  development  of  gasdynamic  lasers,  new  interest 
has  arisen  in  the  study  of  supersonic  flow  diffusers.  The 
steady  state  performance  characteristics  such  as  pressure 
recovery  and  diffuser  length  for  optimum  recovery  of  straight 
duct,  constant  area  and  variable  area  diffusers,  were  most 
recently  investigated  at  our  laboratory  by  Merkli  (1975a,  b, 
1976)  in  a continuous  wind  tunnel.  Starting  times  and  initial 
pressure  ratios  needed  to  start  supersonic  flow  at  the  nozzle 
exit  have  also  been  reported  by  Merkli  and  Abuaf  (1976)  who 
used  a Ludwieg  tube  with  upstream  and  downstream  diaphragm 
locations . 

The  Ludwieg  tube,  or  the  tube  wind  tunnel  (Ludwieg  1955), 
is  an  intermittent  wind  tunnel  where  the  supply  tank  is  replaced 
by  a long  tube.  A converging-diverging  nozzle  generating  super- 
sonic flow  in  a test  section  is  located  between  the  upstream 
and  downstream  sections.  A diaphragm,  separating  the  high 
pressure  side  (p^,  T^)  from  the  low  pressure  one  (p^,  T^), 
before  flow  is  initiated,  can  be  placed  either  upstream  or 
downstream  of  the  test  section.  Ludwieg  used  a downstream 
valve  and  Hottner  (1965)  was  the  first  to  conduct  experiments 
with  an  upstream  diaphragm  location.  Falk  (1963),  Falk  and 
Hertzberg  (1967),  and  Falk  (1968)  presented  an  approximate 
wave  diagram  for  the  starting  process  where  the  nozzle  was 
replaced  by  a zero  length  section,  followed  by  a constant 
area  test  section.  Davis  and  Gwin  (1967)  and  Davis  (1968) 
gave  a similar  wave  diagram  for  the  flow  process  in  a Ludwieg 
tube.  Moreover  they  reported  experiments  with  an  upstream 
diaphragm  location  resulting  in  starting  times  that  were  55 
per  cent  shorter  than  those  obtained  with  a downstream 
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diaphragm  location.  Figure  1 represents  the  wave  diagram  of 
the  Ludwieg  tube  flow  after  diaphragm  rupture  after  Davis 
and  Gwin  (1967)  and  Miller  (1968).  The  general  nature  of  the 
flow  as  indicated  in  this  simplified  diagram  is  basically 
correct,  but  the  wave  process  in  the  throat  region  is  in 
reality  more  complicated  due  to  reflections  of  the  shock 
wave  from  the  converging  walls  of  the  nozzle.  A more  detailed 
wave  analysis  of  the  starting  process  is  in  fact  needed  to 
predict  a more  realistic  wave  diagram  of  the  flow  field. 

The  flow  in  the  Ludwieg  tube  is  initiated  by  the  rupture 
of  the  diaphragm.  In  the  case  of  an  upstream  diaphragm  a 
backward-facing  centered  expansion  fan  propagates  into  the 
supply  tube.  This  fan  is  reflected  off  the  closed  end  of  this 
tube  and  it  returns  to  the  nozzle  throat.  This  process  is 
called  the  first  flow  cycle.  The  expansion  wave  sets  the  gas 
in  motion  and  it  slightly  lowers  the  local  pressure  and 
temperature  (p^,  T^)  ahead  of  the  nozzle.  Since  the  flow  speed 
in  the  nozzle  supply  after  passage  of  the  expansion  fan  is 
low,  p.=  p and  T~  = T , where  the  subscript  o refers  to  the 
nozzle  supply  conditions  in  analogy  to  conventional  super- 
sonic wind  tunnels.  Expressions  relating  all  properties  of 
the  gas  during  the  first  flow  cycle  have  been  derived  by 
Cable  and  Cox  (1963)  and  Davis  and  Gwin  (1967).  Subsequent 
to  the  rupture  of  the  diaphragm,  a shock  wave  and  a contact 
surface  proceed  downstream  through  the  nozzle  and  test  section. 
Additional  waves  follow  the  contact  surface  that  separates 
the  gas  initially  in  the  supply  tube  and  the  downstream  tube, 
to  adjust  the  nozzle  test  section  pressure  to  the  freestream 
pressure  corresponding  to  the  design  Mach  number  expected 
at  the  nozzle  exit  for  a given  nozzle. 

The  total  starting  time,  t^,  at  a fixed  location  is  defined 
as  the  time  that  elapses  between  the  rupture,  of  the  diaphragm 
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and  the  establishment  of  steady  flow,  i.e.  constant  pressure 
at  that  location.  The  starting  time  is  the  sum  of  two  time 
intervals  as  noted  in  Figure  1.  For  one  t ^ gives  the  elapsed 
time  from  the  diaphragm  rupture  until  the  shock  wave  arrives 
at  the  specific  location,  and  t denotes  the  time  between 
the  arrival  of  the  shock  and  the  attainment  of  the  steady 
free  stream  pressure.  This  total  starting  time  definition  is 
dependent  on  location.  At  each  point  of  the  flow  field  a 
different  total  starting  time  will  be  recorded. 

In  the  present  report  we  have  investigated  the  starting 
process  of  the  supersonic  flow  at  the  nozzle  and  diffuser 
exits  respectively  by  means  of  spark  shadowgraph 
pictures  in  which  the  different  phases  of  the  flow  can  be 
seen.  The  results  compare  well  with  the  static  pressure 
measurements  performed  at  the  same  locations.  The  results 
of  the  flow  starting  times  and  corresponding  initial  pressure 
ratios  needed  to  start  the  flow  at  the  nozzle  exit,  are 
extended  to  a wider  range  of  initial  upstream  pressures  and 
Reynolds  numbers  than  the  results  previously  reported  (Merkli 
and  Abuaf  1976).  Moreover  the  effects  of  boundary  layer 
suction  at  the  nozzle  exit  on  starting  time  was  determined. 

II.  EXPERIMENTAL  METHODS 

The  experimental  investigation  was  conducted  in  the  Yale 

Ludwieg  tube  which  was  described  in  detail  in  our  last  report 

(Merkli  and  Abuaf  1976).  The  upstream  tube  (5.3  inches  ID, 

22  feet  long)  and  the  downstream  lube  (3.76  inches  ID,  25  feet 
3 

long,~  2 ft  ) sections  respectively  have  9CE  bends  to  accommodate 

the  total  tube  length  within  the  confines  of  the  laboratory 

(Fig.  2a).  The  downstream  section  was  connected  to  a dump 
3 

tank  (17  ft')  to  reduce  the  effect  of  shock  wave  reflections. 
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The  24  inch  long  test  section,  (2x5)  in  , is  connected  to 
both  the  upstream  and  downstream  sections  by  two  transition 
pieces.  For  the  reported  experiments  the  diaphragm  ( 3 mil 
Mylar  sheets)  was  located  at  the  upstream  end  of  the  test 
section  ( 1.81  inches  from  the  nozzle  entrance).  The  dia- 
phragms were  scratched  along  diagonal  lines  with  a razor 
blade  and  they  were  ruptured  by  an  X-shaped  cutter.  This 
preparation  prevents  diaphragm  pieces  from  flying  through 
the  test  section. 

A 12  inch  long  converging-diverging  nozzle  is  placed 
in  the  test  section,  leaving  another  12  inches  of  test  sec- 
tion for  the  installation  of  various  diffusers.  The  two- 
dimensional  nozzle  contour  was  based  on  method  of  charac- 
teristics calculations  for  a uniform  M=3  flow  at  the  nozzle 

exit  without  a boundary  layer  correction  (Merkli  1976).  The 

. 2 

throat,  (0.472x2)  m , is  located  at  4 inches  from  the 

2 

nozzle  inlet  and  the  nozzle  exit  area  is  (2x2)  in  . Metal 
side  walls  with  holes  at  1 inch  intervals  along  the  flow 
axis  allow  pressure  measurements  to  be  performed.  For  these 
static  pressure  measurements  Kistler  Model  606L  quartz 
pressure  transducers  (rise  time  3 ps)  are  used  with  a Kistler 
Model  504  charge  amplifier. 

After  placing  the  diaphragm, in  a typical  experiment  both 
the  upstream  and  downstream  sections  are  pumped  down  to  a few 
torr.  The  upstream  section  is  next  filled  with  dry  air  (dew 
point  better  than  213  K)  to  a given  initial  upstream  pressure, 
p^ , which  is  measured  either  by  a Wallace  & Tiernan  Model 
FA145  pressure  gauge  (0-120  inches  of  mercury,  ± 0.2  inches 
Hg)  or  by  a Wallace  & Tiernan  Type  FA187  mercury  manometer, 
capable  of  measuring  an  accuracy  of  0.1  torr  in  a range  of 
0-800  torr.  The  initial  downstream  pressure,  p^ , is  next  set 
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to  the  required  value  by  bleeding  air  from  the  atmosphere 
and  it  is  measured  by  Wallace  S Tiernan  Model  FA160  gauges 
(0-100  torr  or  0-400  torr) . The  breaking  of  the  diaphragm 
initiates  the  flow  through  the  test  section.  The  output  of 
a pressure  transducer  located  at  the  nozzle  entrance  (x=-ll 
inches)  triggers  a four  channel  storage  oscilloscope 
(Tektronix  Model  5103N).  Typical  oscillograms  of  the  pressure 
variation  as  a function  of  time  are  given  in  Fig.  2b. 

Two  plexiglass  windows  allowed  shadowgraph  and  schlieren 
observations  of  the  entire  nozzle-diffuser  flow  area.  For  the 
shadowgraph  pictures  of  the  flow  field  a spark  source  was 
used  with  a collimating  lens  (4  inches  in  diameter)  as  noted 
in  Fig.  2a.  The  output  of  a pressure  transducer  located  at  the 
throat*  and  sensing  the  arrival  of  the  shock  wave,  triggered 
a Hewlett  Packard  5233L  electronic  counter  and  a General 
Radio  Corp . Type  1392-A  time  delay  generator.  After  a preset 
time  delay,  an  output  from  the  time  delay  generator  triggered 
the  short  duration  high  voltage  (1  ps  , 7500  v)  spark  source. 

The  output  of  a light  sensitive  phototransistor  was  used  to 
stop  the  electronic  counter  at  the  occurrence  of  the  spark. 

Thus  the  time  period  between  the  arrival  of  the  starting  shock 
at  the  nozzle  throat  and  the  recording  of  the  photograph  can 
be  directly  read  from  the  electronic  counter  within  an 
accuracy  of  1 ps. 

III.  EXPERIMENTAL  RESULTS 

1.  STARTING  PROCESS 

* The  pressure  transducer  to  trigger  the  counter  and  the  time 
delay  generator  was  located  at  the  throat  because  at  that  loca- 
tion a well  defined  shock  wave  and  pressure  pulse  appear.  This 
scheme  eliminates  the  uncertainties  involved  in  the  non-repro- 
ducible breaking  of  the  diaphragm. 
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When  the  upstream  diaphragm  of  a Ludwieg  tube  is  ruptured 
a wave  pattern  as  presented  in  Figure  1 develops  in  the  nozzle 
followed  by  the  straight  duct  diffuser.  The  present  experimental 
results  generally  agree  with  this  picture. 

A.  Shadowgraph  pictures  of  the  starting  process 

Subsequent  to  the  rupture  of  the  diaphragm  a shock  wave 
travels  along  the  nozzle-diffuser  duct  (Fig.  1).  This  fact  is 
observed  in  the  typical  shadowgraphs  presented  in  Fig.  3 taken 
at  different  locations  and  times.  An  upstream  supply  pressure, 

Pq  = 966  torr , an  initial  downstream  pressure,  p-^  = 150  torr,  were 
used  with  the  nozzle  followed  by  a 6 inches  long  diffuser  of 
constant  cross  section.  The  photographs  are  given  at  full  scale. 
The  vertical  solid  black  lines*  are  strings  attached  at  reference 
points  along  the  test  section  in  order  to  determine  the  location 
of  the  shock  wave.  The  first  picture  shows  the  pressure  trans- 
ducer located  at  the  throat  and  the  last  picture  in  the  sequence 
depicts  the  shock  leaving  the  end  of  the  diffuser.  The  bound- 
ary layer  estimated  to  be  less  than  0.09  inches  at  the  nozzle 
exit  for  the  present  experimental  conditions  is  not  visible 
in  the  pictures.  The  x-t  diagram  for  the  shock  wave,  for  the 
initial  pressure  ratio,  p4/p^=6.67  is  presented  in  Fig.  4 as 
obtained  from  similar  shadowgraph  pictures.  In  the  same  figure, 
a theoretical  curve  is  plotted  for  comparison.  First  an 
initial  shock  Mach  number  was  calculated  for  the  given  pres- 
sure ratio,  p|+/p^  = 6.67,  by  using  the  following  relationship 
(Glass  and  Hall  1959) 

p4/p1  = 1/6  [7M^-l].[l-l/6(Ms-l/Ms)]'7  , 


* All  the  extra  lines  present  in  the  picture  are  caused  by 
scratches  on  the  plexiglass  windows. 
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valid  for  the  case  of  air/air  and  a constant  area  shock 
tube.  The  variation  of  the  shock  Mach  number  with  the  cross 
sectional  area  as  it  travels  through  the  converging-diver- 
ging nozzle  was  taken  into  account  by  the  following  expres- 
sion (Chisnell  1957,  Whitham  19 58 ) 

dA  2 M dM 

_ _ s s 

A (M2-l)  K(M  ) 

s s 

An  expression  for  the  function  KCM^)  is  given  in  the  above 
references.  The  x-t  diagram  for  the  shock  wave  is  drawn 
knowing  the  speed  of  sound  for  T^=296  K and  the  variation  of 
the  shock  Mach  number,  M , along  the  nozzle.  The  experimental 
results  recorded  downstream  of  the  throat  agree  well  with  the 
theoretical  predictions  with  a maximum  deviation  of  12  %. 

Shadowgraph  pictures  were  also  taken  in  the  vicinity  of 
the  6 inch  diffuser  exit  for  various  initial  pressure  ratios, 
p^/p^,  recording  the  location  of  the  shock  wave  and  the  time 
it  took  for  the  shock  wave  to  travel  from  the  nozzle  throat 
to  the  specific  location.  This  provided  an  average  experi- 
mental shock  wave  Mach  number,  M . The  results  are  depicted 
in  Fig.  5 where  a solid  line  gives  the  shock  Mach  number  as 
predicted  by  the  constant  area  shock  tube  theory  (Glass  and 
Hall  1959).  Moreover  a dashed  line  gives  an  average  shock 
Mach  number  calculated  by  taking  into  account  the  area  change 
along  the  nozzle  (Whitham  1958).  As  observed  in  Fig.  6,  the 
experimental  results  show  scatter,  however  they  are  still 
within  ± 10  % of  the  latter  predictions.  The  experimental 
average  shock  Mach  numbers  are  lower  than  the  predicted  ones 
for  the  low  pressure  ratios.  Conversely  they  are  higher  for 
the  higher  ones.  The  speed  of  the  diaphragm  rupture  depends 
strongly  on  the  initial  pressure  ratio.  Since  the  nozzle 
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entrance  is  very  close  to  the  diaphragm  location  (at  1.81 
inches),  the  effect  of  the  diaphragm  rupture  will  therefore 
affect  the  experimental  results. 

According  to  Fig.  1,  the  upstream  facing  wave* 
the  interface  will  adjust  the  pressure  level  at  a given 
location  to  the  steady  free  stream  pressure  which  corresponds 
to  the  local  Mach  number.  Fig.  6 presents  four  typical 
shadowgraphs  taken  with  p1+/p-^  = 6.67  at  different  times  for 
the  same  test  section  configuration.  The  pictures  show  the 
upstream  facing  wave  which  is  formed  by  the  effects  on  the 
flow  by  the  interaction  of  the  initial  pressure  ratio  and 
the  local  steady  free  stream  pressure.  The  wave  as  expected 
becomes  stronger  and  thus  more  visible  downstream  of  the 
nozzle  throat.  This  wave  moves  then  downstream  and  stops 
at  a given  location  along  the  test  section  which  is  mainly 
determined  by  the  initial  pressure  ratio,  p^/p^.  The  x-t 
diagram  for  this  wave  as  obtained  from  similar  shadowgraph 
pictures  is  shown  in  Fig.  7.  The  time  measured  here,  t ^ , 
corresponds  to  the  time  interval  between  the  arrival  of  the 
shock  at  the  nozzle  throat  and  the  occurrence  of  the  spark, 
i.e.  the  instant  at  which  the  picture  is  taken.  The  x-t 
diagram  for  the  starting  shock  is  also  plotted  on  the  same 
figure  for  comparison.  The  results  agree  qualitatively  with 
the  wave  pattern  presented  in  Fig.  1.  At  this  stage,  these 
results  for  the  upstream  facing  wave  can  not  be  compared  with  a 
theoretical  curve,  since  no  such  calculations  are  available. 

A series  of  shadowgraphs  taken  at  the  exit  of  the  nozzle 
as  a function  of  time  are  presented  in  Fig.  8.  The  vertical 
black  line  is  again  a string  attached  at  0.2  inches  down- 
stream of  the  nozzle  exit.  First  the  shock  is  observed, 
then  at  3 ms  some  disturbances  pass  through,  and  finally  at 

* By  an  upstream  facing  wave,  we  mean  a Q shock  moving  down- 
stream and  being,  accompanied  by  a local  pressure  decrease 
after  its  passage. 
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4 to  7 ms  the  upstream  facing  wave  has  reached  the  nozzle 
exit  and  it  is  observed  to  oscillate  back  and  forth  about 
the  flow  direction  at  that  location.  Steady  supersonic 
flow  with  free  stream  conditions  follows  this  upstream 
facing  wave.  The  wave  system , following  the  initial  shock 
wave  and  ending  with  the  upstream  facing  wave,  appears  as 
oblique  shocks  at  the  nozzle  exit,  followed  by  a series  of 
shock  waves  within  the  diffuser.  Deceleration  of  the  flow, 
adverse  pressure  gradient  effects,  boundary  layer  separation 
and  turbulent  mixing  are  also  observed. 

When  the  upstream  facing  wave  reaches  a given  location 
in  the  nozzle-diffuser  assembly  oscillating  around  that 
point,  steady  supersonic  flow  ceases  at  that  location  and 
the  flow  breaks  down  (Merkli  and  Abuaf  1976).  When  the 
upstream  facing  wave  stays  downstream  of  the  nozzle  exit 
for  p^/p,  > 6.67  steady  supersonic  flow  exists  throughout 
the  nozzle  and  it  is  followed  by  a slowing  down  and 
pressure  increase  in  the  diffuser  through  a series  of  shock 
waves.  The  corresponding  pressure  rise  yields  the  pressure 
recovery  and  the  length  of  the  shock  system  in  the 
diffuser  is  called  the  pressure  recovery  zone  (Merkli 
1975a,  1976).  Varying  the  initial  pressure  ratio,  p1+/p-^, 
for  a 6 inch  diffuser,  the  location  where  the  upstream 
facing  wave  stops,  and  where  the  supersonic  flow  breaks 
down  in  the  nozzle-diffuser  assembly  can  be  determined 
(see  Fig.  9).  For  p^/p-^  > 10.1,  supersonic  flow  is 

observed  throughout  the  entire  test  section  with  the 
familiar  over-  and  under-expanded  jet  configurations  at 
the  diffuser  exit.  For  10.1  > pl(/p^  > 6.67,  the 

disturbance  moves  from  the  diffuser  end  upstream  towards 
the  nozzle  end.  For  p1(/p^  < 6.6  7 the  shock  wave  system 
moves  further  upstream  into  the  nozzle.  Since  the  upstream 
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facing  wave  is  observed  to  oscillate  around  a given  location, 
the  limiting  line  as  presented  in  Fig.  9 for  the  break  down 
of  the  flow  will  in  reality  be  given  by  a broader  range  of 
conditions . 

B.  Starting  times 

The  total  starting  time  at  a given  location  is  defined  as 

the  time  period  required  to  establish  steady  flow  at  that 

location  after  the  rupture  of  the  diaphragm.  As  discussed, the 

breaking  time  of  the  diaphragm  is  erratic  and  thus  a source 

of  error.  To  eliminate  this  uncertainty  in  shock  tube  research, 

it  is  customary  to  assume  that  this  time  is  negligibly  small 

(Fig.  1).  The  previous  section  was  concerned  with  the  ti/ne 

period  between  the  rupture  of  the  diaphragm  and  the  arrival 

of  the  shock  wave  at  a given  location,  t , . Its  variation 

’si 

was  investigated  as  a function  of  the  experimental  parameters 

and  an  approximate  method  for  its  prediction  was  presented. 

This  section  will  deal  with  the  time  required  to  establish 

steady  flow  after  the  shock  wave  has  arrived  at  that  location, 

i.e.  t^ 2 • The  results  are  similar  to  those  presented  in  the 

previous  report  by  Merkli  and  Abuaf  (1976)  except  that  they 

are  extended  to  a wider  range  of  supply  pressures  or 

Reynolds  numbers  based  on  the  properties  calculated  for  the 

free  stream  conditions  at  the  nozzle  exit. 

Fig.  10  depicts  the  second  period  of  the  flow  starting 

time,  t 0,  at  the  nozzle  exit  without  a diffuser  as  a function 
s 2 

of  the  pressure  ratio,  P^/po  for  various  nozzle  supply 
pressures,  p^  or  Reynolds  numbers,  Re„  *.  The  wide  scatter 
of  the  experimental  data  is  not  caused  by  errors  in  measurement 

*Here  the  length  entering  the  Reynolds  number  is  defined  as 
the  hydraulic  diameter  of  the  nozzle  exit,  D=2  inches. 
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but  rather  it  is  due  to  the  erratic  breaking  of  the  Mylar 

diaphragms.  This  scatter  diminishes  for  the  data  recorded 

at  the  diffuser  exit  for  longer  diffusers.  The  starting 

times,  t_  2 5 increase  linearly  with  the  pressure  ratio,  p^/pQ, 

and  it  is  independent  of  the  upstream  supply  pressure  or 

c.  5 

Reynolds  number  for  2.8x10  < Re^  < 7.2x10  . The  vertical 

line  at  the  right  hand  side  of  the  figure  corresponds  to  the 
limiting  pressure  ratio,  p^/pQ,  for  which  the  steady  super- 
sonic flow  at  the  nozzle  exit  was  disturbed.  Figs.  11  and 

12  present  the  flow  starting  times,  t at  the  nozzle  and 
diffuser  exits  respectively  for  a 3 inch  long  diffuser.  Figs. 

13  and  14  depict  similar  results  for  a 6 inch  diffuser.  Here 
too  the  starting  times,  t £>  increase  linearly  with  the 
pressure  ratio,  p./p  , and  they  are  independent  of  the  supply 

1 ° 5 

pressure,  p , or  Reynolds  number,  Re^,  for  3.8x10  < Re^  < 

1.7x10^.  Figs.  15,  16  and  17,  18  finally  present  similar 
data  for  a 9 and  a 12  inch  long  diffuser  respectively.  Fitting 
a straight  line  through  the  experimental  data  of  each  case 
by  means  of  a least  square  analysis  provides  the  results 
tabulated  in  Tables  I and  II  respectively.  The  straight  lines 
obtained  from  these  expressions  are  the  solid  lines  shown 
in  Figs.  10  to  18. 

A cross  plot  of  the  starting  times,  tt^  > for  the  flow  at 
the  nozzle  and  diffuser  exits  respectively  as  a function  of 
diffuser  length  for  a constant  value  of  the  initial  pressure 
ratio,  p^/pQ,  gives  the  curves  presented  in  Figs.  19  and  20. 
For  low  values  of  the  pressure  ratio,  p^/p  < 0.05  the 
starting,  time  for  the  flow  at  the  nozzle  exit  is  nearly 
independent  of  the  diffuser  length.  Tor  higher  values  of 
the  pressure  ratio,  p^/p  > 0.075,  the  flow  starting  times 
at  the  nozzle  exit,  t increase  with  the  diffuser  length. 
Fig.  20  depicting  the  variation  of  the  flow  starting  times, 
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ts2»  at  the  diffuser  exit  shows  that  for  a given  value  of 
the  pressure  ratio,  p^/po,  the  flow  starting  times  increase 
linearly  with  the  diffuser  length.  The  variation  of  the 
dimensionless  starting  times,  t „ = t with  the  axial 

distance,  x/D,  are  plotted  in  Figs.  21  to  24  for  various 
length  diffusers,  L/D=  1.5,  3,  4.5,  and  6,  and  for  different 
values  of  the  pressure  ratio,  p-^/p^  The  characteristic  flow 
time  here  defined  by  T^=l/a*,  is  related  to  the  geometry  of 
the  nozzle  throat  and  it  is  a measure  of  the  acceleration 
of  the  flow  or  its  cooling  rate  (Wegener  and  Cagliostro  1972). 
The  characteristic  length,  Z , =/  R*h* , where  h*  is  the  height 
of  the  nozzle  at  the  throat  (h*=0.472  inches),  and  R*  is  the 
radius  of  curvature  of  the  nozzle  profile  at  the  throat  (R*= 
7.64  inches).  The  speed  of  sound  at  the  sonic  throat  is  given 
by  a*.  At  locations  where  x/D  < -3  the  pressures  recorded  by 
the  pressure  transducers  show  an  increase  in  the  pressure 
level  above  the  initial  value  with  the  arrival  of  the  shock. 
The  pressure  overshoots  its  steady  state  level,  reaches  a 
maximum,  and  then  drops  down  and  levels  off  at  its  steady 
state  running  value  (see  Fig.  2b).  The  starting  time,  t,^,  at 
those  points  was  taken  as  the  time  period  between  the  arrival 
of  the  shock  and  the  first  time  the  pressure  level  reached 
the  steady  state  running  value.  In  this  region  the  flow 
starting  times  are  seen  to  be  independent  of  the  pressure 
ratio,  p^/pQ,  and  °T  "the  diffuser  length.  In  the  region 
-3  < x/D  < -2  a sharp  increase  in  the  starting  times, 
is  observed.  This  is  the  same  location  where  the  upstream 
facing  wave  was  observed  in  Figs.  6 and  7.  This  region  is 
seen  to  move  upstream  closer  to  the  nozzle  throat  for  the 
longer  diffusers,  i.e.  from  x/D=-2  for  L/D=1.5  to  x/D=-3 
for  L/D=6. 

C.  Minimum  initial  pressure  ratios  needed  to  start  the 


flow . 
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The  initial  pressure  ratios,  p^p.^  , sufficient  to  start 
the  supersonic  flow  at  the  nozzle  and  diffuser  exits  are 
plotted  in  Fig.  25  as  a function  of  the  diffuser  length, 

0 <L/D  <6.  The  pressure  ratios,  pQ/p^,  that  start  the  flow 
at  the  nozzle  exit,  show  a minimum  (p  /p^=5.2)  once  a 
diffuser  length  of  L/D=4.5  is  attained.  In  turn  the  optimum 
pressure  ratio,  p^/p^ , sufficient  to  start  the  flow  at  the 
nozzle  exit  without  a diffuser  is  equal  to  p /p^=10.3  (±  10%). 

It  is  interesting  to  note  that  no  change  of  the  needed 
starting  pressure  ratio  for  the  nozzle  exit  without  a dif- 
fuser occurs  when  diffusers  of  varying  length  are  added 
whose  exit  flow  is  supersonic.  In  fact  the  diffuser  acts 
here  like  an  extended  parallel  wall  test  section.  The  pre- 
sent results  agree  with  the  conclusions  of  our  last  report 
(Merkli  and  Abuaf  1976  ) and  the',  are  independent  of  the 

Reynolds  number  in  the  restricted  range  2.8x10^  < Ren  < 

6 ^ 

1.7x10  which  was  used  in  this  investigation. 


2.  PRESSURE  RECOVERY  DURING  THE  FIRST  FLOW  CYCLE  OF 
THE  LUDWIEG  TUBE 

The  static  pressure  levels  at  the  diffuser  exit,  p^ , for 
diffusers  of  different  lengths  were  also  recorded  during  the 
first  flow  cycle  of  the  Ludwieg  tube  while  operating  at  the 
limiting  pressure  ratios,  p^/p^.  Under  these  conditions 
the  supersonic  flow  at  the  nozzle  exit  was  just  starting 
to  be  disturbed  (see  Fig.  25).  The  pressure  ratio,  pQ/p^, 
which  should  correspond  to  the  steady  state  pressure  recovery 
of  the  diffuser,  is  plotted  in  Fig.  26  for  various  diffuser 
lengths.  The  pressure  levels  at  the  diffuser  exits  showed 
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violent  fluctuations  with  lime  during  the  first  flow  cycle 
(up  to  ± 20-30  % of  p^)  leading  to  the  uncertainty  indicated 
by  the  error  bars  (Fig.  26).  The  square  symbol  presents  the 
point  obtained  by  Merkli  (1976)  with  a steady  continuous 
wind  tunnel  for  a similar  nozzle  and  for  the  boundary  layer 
displacement  parameter  <5S'C/D  of  about  0.04.  The  current  results 
agree  with  those  of  the  continuous  tunnel.  The  pressure 
fluctuations  at  the  diffuser  exit  during  the  steady  operation 
in  the  first  cycle  at  the  pressure  ratio  when  the  flow  at 
the  nozzle  exit  is  breaking  down  may  be  caused  by  the 
following.  For  one  , the  operation  period  of  the  Ludwieg  tube 
(first  cycle)  is  not  as  steady  as  predicted  from  theory  since 
variations  of  the  downstream  pressure  by  reflection  of  waves 
may  cause  pressure  fluctuations.  It  is  well  known  that 
close  to  optimum  performance  diffuser  flows  are  very  sensitive 
to  small  disturbances.  Secondly  the  pressure  variations  may 
be  due  to  adverse  pressure  gradient  effects.  Boundary  layer 
separation  and  turbulent  mixing  may  be  carried  downstream 
to  the  transducer  location  at  the  diffuser  exit.  The  continuous 
wind  tunnel  experiments  (Merkli  1976)  were  carried  out  with 
a Statham  gauge  which  is  insensitive  to  fast  pressure 
fluctuations  and  thus  the  fluctuations  may  have  been  overlooked. 

3.  FLOW  STARTING  WITH  BOUNDARY  LAYER  SUCTION  AT  THE  NOZZLE 
EXIT. 

To  investigate  the  effect  of  boundary  layer  suction  at 
the  nozzle  exit  on  the  starting  of  the  flow,  a special  3 inch 

long  diffuser  block  was  designed  (see  Fig.  2a  lower  left).  A 

. . 2 
suction  slit  consisting  of  a (3/16x1.5)  in  opening  located 

at  the  nozzle  exit  was  added.  The  dump  tank  was  disconnected 

from  the  downstream  section  of  the  Ludwieg  tube  and  the  90° 

bend-T  piece  was  replaced  by  a simple  90°  bend.  The  new 


diffuser  blocks  were  connected  by  copper  tubing  (1  inch  ID, 

7 feet  long)  to  the  dump  tank  through  a second  diaphragm 
and  a pin  mechanism.  To  conduct  an  experiment,  the  initial 
upstream  and  downstream  pressures,  p^  and  p^  , respectively 
were  set  to  the  required  value  following  the  procedure 
described  in  Chapter  II.  The  large  dump  tank  was  evacuated. 
First,  the  diaphragm  separating  the  dump  tank  from  the  copper 
tubing  leading  to  the  diffuser  block  was  ruptured  in  order 
to  intiate  suction  through  the  slit.  Next  the  upstream  dia- 
phragm of  the  Ludwieg  tube  was  ruptured  manually  starting  the 
flow  through  the  test  section.  The  flow  was  assumed  to  be  choked 
at  the  diffuser  suction  slit.  For  a given  gas  and  a constant 
supply  temperature,  the  mass  flow  rate,  m,  is  proportional 
to  the  throat  area  and  the  supply  pressure  for  choked  flow. 

The  pressure  at  the  exit  of  the  nozzle  during  the  first  flow 
cycle  is  about  35  torr,  (corresponding  to  a supply  pressure, 

PQ , of  966  torr),  and  the  area  of  the  two  slits  is  0.56  in  . 

The  suction  mass  flow  rate  is  thus  around  2 % of  the  total 
mass  flow  rate  through  the  nozzle  in  the  first  cycle.  Starting 
times  for  the  flow,  t „ , at  the  nozzle  and  diffuser  exits 
with  a 3 inch  diffuser  and  suction  are  presented  in  the  last 
two  Figs.  27  and  28  respectively  as  a function  of  the  initial 
pressure  ratio,  p^/p^  It  is  clear  that  the  starting  times 
at  both  locations  are  reduced  by  16  to  33  % by  suction  for 
the  high  and  low  values  of  the  initial  pressure  ratio,  Pj/P0, 
respectively,  when  we  compare  them  to  the  flow  starting 
times  without  suction  at  the  same  location.  For  the  3 inch 
diffuser  the  boundary  layer  suction  did  not  seem  to  affect 
the  limiting  initial  pressure  ratio,  p /pQ,  which  caused 
the  steady  supersonic  flow  at  the  nozzle  exit  to  be  disturbed. 

One  problem  encountered  during  these  runs  was  the  sequence 
of  the  timing  between  the  initiation  of  the  suction  at  the 
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nozzle  exit  and  the  starting  of  the  flow  in  the  test  section. 
Two  diaphragm  breaking  mechanisms  were  needed  whose  opening 
periods  are  reproducible  with  an  accuracy  of  1 ms.  Thus  the 
expansion  fan  of  the  suction  and  the  initial  shock  wave  of 
the  main  flow  arrive  at  the  nozzle  exit  at  the  same  time. 

The  repeatability  could  not  be  improved  to  better  than  10  ms 
unless  a new  design  would  have  been  introduced.  Thus  in  some 
experiments  the  set  initial  downstream  pressure  was  observed 
to  change  even  before  the  Ludwieg  tube  flow  started. 

IV.  CONCLUSIONS 

1.  The  starting  time  period  from  the  rupture  of  the 
diaphragm  until  the  arrival  of  the  shock  at  a given  location, 
t , can  be  predicted  with  an  accuracy  of  ± 10%  by  predicting 
the  shock  strength  from  the  initial  pressure  ratio,  p^/p  , 
and  by  taking  into  account  the  effect  of  area  change  on  the 
shock  strength  as  it  propagates  through  the  test  section. 

2.  The  starting  time  period  from  the  arrival  of  the  shock 

until  steady  flow  is  established,  t^2>  at  the  nozzle  and 

diffuser  exits;  i)  increases  linearly  with  the  initial  pressure 

ratio,  p./p  , for  a given  diffuser,  ii)  is  independent  of 
a o S 6 

the  Reynolds  number  in  the  range  2.8x10"  _o  1.7x10°,  and  iii) 

for  the  initial  pressure  ratio  p-^/p  < 0.05,  the  flow 

starting  time  at  the  nozzle  exit,  t 2,  is  independent  of  the 

diffuser  length,  but  increases  linearly  with  diffuser  length 

for  higher  pressure  ratios. 

3.  The  initial  pressure  ratio,  pQ/p^,  to  start  the  super- 
sonic flow  at  the  nozzle  exit  shows  a minimum  (p  /n,=5.2) 
for  a diffuser  length  of  L/D=4.5,  and  it  is  independent  of 
the  Reynolds  number,  Re^,  in  the  range  quoted. 

4.  The  pressure  ratio,  PQ/p(j  , during  the  first  flow  cycle 
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of  the  Ludwieg  tube,  corresponding  to  the  pressure  recovery 

through  the  diffuser,  shows  a minimum  value  (p  /p  =4.45)  for 

^o  1 d 

a diffuser  length  of  L/D=4.5.  These  results  are  in  agreement 
with  those  obtained  by  Merkli  (1976),  L/D  ~ 6 and  po/p^=4.2, 

in  a continuous  wind  tunnel  under  similar  experimental 
conditions . 

5.  For  a 3 inch  diffuser  with  boundary  layer  suction  to 

remove  2%  of  the  total  mass  flow  rate  through  the  nozzle 

applied  through  a slit  at  the  nozzle  exit  appreciably  reduced 

the  starting  times  for  the  flow,  t „ , at  the  nozzle  and 

s 2 

diffuser  exits. 


The  supersonic  diffuser  performance  characteristics 
observed  can  be  summarized  as  follows.  First  the  simplest 
geometry  of  constant  area  emerged  as  quite  attractive,  efficient, 
and  easy  to  construct.  The  supersonic  flow  starting  times 
at  the  nozzle  exit  are  shorter  for  an  upstream  diaphragm 
location  than  a downstream  one.  These  starting  times  can 
still  be  reduced  by  boundary  layer  suction  at  the  nozzle 
exit.  The  flow  starting  process  can  in  principle  be 
predicted  from  wave  diagrams  existing  in  the  literature. 

One  disadvantage  of  this  configuration  is  the  diaphragm 
material  flying  through  the  test  section  which  can  disturb 
the  starting  flow  field,  but  this  can  be  eliminated  by 
proper  design  considerations.  The  initial  pressure  ratios 
needed  to  start  the  supersonic  flow  throughout  the  nozzle 
are  not  affected  by  the  location  of  the  diaphragm.  The 
optimum  pressure  recovery  during  steady  operation  and  the 
diffuser  length  required  to  obtain  this  recovery  can  be 
predicted  once  the  Mach  numbei'’  and  the  boundary  layer 
displacement  parameter  are  known  at  the  exit  ot  the  nozzle. 

Although  we  are  far  from  a complete  undo)  . 1 Hiding  of 
supersonic  diffuser  flows,  the  experimental  findings  reported 
are  providing,  some  basic  points  to  consider  in  the  choice  and 
design  of  supersonic  diffusers,  for  some  specific  applications 
such  as  gasdynamic  lasers. 
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TABLE  I 

VARIATION  OF  THE  STARTING  TIME,  t^tms),  WITH  THE  INITIAL 
PRESSURE  RATIO,  p1/po,  FOR  THE  NOZZLE  EXIT,  x = 0 inches, 

FOR  VARIOUS  DIFFUSERS.  UPSTREAM  DIAPHRAGM  LOCATION. 

NOZZLE  t _ = 31.7  (p./p  ) + 0.06  0<p./p  <0.1 

s2  f 1 ro  1 1 1 o 

NOZZLE  + 3 inch  DIFFUSER  t 0 = 17.87  (p,/p  ) + 0.63  0<p./p  <0.13 

s2  ^1  1 o r 1 1 o 

NOZZLE  + 6 inch  DIFFUSER  t „ = 23.3  (p./p  ) + 0.50  0<p. /p  <0.15 

s 2 1 1 1 o 1 o 

NOZZLE  + 9 inch  DIFFUSER  tg2  = 28  (Pj/Pq)  + 0.32  0<p1/pQ<0.18 

NOZZLE  + 12  inch  DIFFUSER  t = 38.16  (p./p  ) - 0.25  0<p./p  <0.20 

s2  ^1*0  1 o 


TABLE  II 

VARIATION  OF  THE  STARTING  TIME,  tt,?(ms),  WITH  THE  INITIAL 

PRESSURE  RATIO,  P;L/po,  FOR  THE  DIFFUSER  EXIT  TOR  VARIOUS 

DIFFUSERS.  UPSTREAM  DIAPHRAGM  LOCATION. 

NOZZLE  t „ = 31.7  (p./p  ) + 0.06  0<p./p  <0.1 

(at  x = 0 inches) 

NOZZLE  + 3 inch  DIFFUSER  t „ = 28.76  (p  /p  ) + 0.47  0<p  /p  <0.1 
(at  x = 3 inches) 

NOZZLE  + 6 inch  DIFFUSER  t = 38.66  (p./p  ) + 0.26  0<p,/p  <0.1 

/■  _ _ _ r u \ S l 


NOZZLE  + 9 inch  DIFFUSER  t „ 
(at  x = 9 inches) 


52.97  (p  /p  ) + 0.23  0<p1/po<0.1 


NOZZLE  + 12  inch  DIFFUSER  t „ - 54.01  (p./p  ) + 0.51  0<p./p  <0.1 
(at  x = 12  inches)  S?  1 ° 1 ° 


J 
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Figure  2a.  The  experimental  set-up. 
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